LAL 98-67 
31st August 1998 



Extraction of a s and constraint on the 
Higgs mass from Electroweak fits at 
the Z resonance 



Edwige TOURNEFIERQ 
Laboratoire de l'Accelerateur Lineaire 
IN2P3-CNRS et Universite de Paris-Sud, 
BP34, F-91898 Orsay Cedex 



Abstract 
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1 Introduction 



From 1990 to October 1995 the LEP e + e~ storage ring was operated at 
center-of-mass energies close to the Z mass (called LEP1 program). LEP1 
data have been collected and analysed by the 4 LEP experiments ALEPH(l|], 
DELPHI [|, L3@ and OPALg], the electroweak results and their combina- 
tion are given in Ref. [||] and are still preliminary Since 1992 the SLD detector 
is taking data on the SLC e + e~ storage ring operating at center-of-mass en- 
ergies also close to Mz with polarised electron beam. The preliminary results 
from SLD include 1992 to 1997 data §. 

The cross sections and the asymmetries of the reactions e + e~ — > ff (7) mea- 
sured at LEP and SLC are sensitive, through radiative corrections, to the 
following Standard Model parameters: the strong coupling constant a s , the 
Higgs mass Mh and the top mass M t . The value of a s is mainly determined 
by the LEP total cross section measurements, while the asymmetries mea- 
sured at LEP and SLC are most sensitive to Mh- We will concentrate on the 
determination of a s and therefore on the determination of the Z lineshape 
parameters at LEP. 

2 The Z lineshape and the fitting procedure 

At LEP, an integrated luminosity of 110 pb _1 per experiment has been accu- 
mulated at the Z peak and about 40 pb _1 off peak, mostly at Mz ±1.8 GeV. 
About 4 million hadronic and 0.5 million leptonic Z decays have been col- 
lected by each of the four LEP experiments. This large sample allows a 
precise determination of the Z boson properties: the Z mass Mz, the Z width 
Tz, the total hadronic cross section at the pole <r° ad and the ratio of hadronic 
to leptonic pole cross sections = cr^/of = T hSud /Ti. 
These parameters have the advantage of being almost uncorrelated: 

• Mz is determined by the position of the peak and therefore depends on 
the absolute energy scale. Its uncertainty is dominated by the uncer- 
tainty on the LEP energy which is about 1.5 MeVf|. 

• Tz is determined by the width of the Z resonance and therefore by 
peak and off peak relative cross section measurements. The uncertainty 
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on Tz comes mainly from uncorrelated errors on the off peak energy 
measurement (~1.5 MeV0]) and from off peak statistics. 

• ^had i s determined by the height of the resonance and is derived from 
the measurement of the hadronic cross section. The statistical and 
systematic uncertainties of the selections are typically of the order of 
0.5 x and 0.8 x 10~ 3 respectively, for each experiment 0. The 
error on <r° ad is dominated by the theoretical uncertainty on the lumi- 
nosity which is 1.1 x 10~ 3 . Note that a recent study on the theoretical 
precision of the LEP luminosity |7j] will lead to a reduction of this error 
to 0.6 x 10" 3 . 

• R,£ is determined by the measurement of the leptonic and hadronic cross 
sections. Since is a ratio, the uncertainty arising from the luminosity 
cancels and thus the uncertainty on R^ comes only from the statistical 
and systematic errors in the e + e~ — > ff event selections which will be 
discussed in Section |3|. 



Figure 1: Lowest order diagrams contributing to the e + e — > ff cross section. 

The lowest order diagrams involved in the process e + e~ — > ff are shown in 
Figure [l]. The cross section is the sum of the Z exchange, the 7 exchange 
and their interference 

a = a z + o-y + (T illt (1) 

The cross section due to the Z exchange is parametrised with a Breit-Wigner 
in a model independent way using Mz, Tz, and the Z partial widths Tf 



u sr z -,i i2vrr e r f 

ffz = a < (s - m + (ir|7M|? Wlth °> = ^T (2) 
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Parameter 


LEP average 


M z 


91.1867 ±0.0020 


r z 


2.4948 ± 0.0025 


°had 


41.486 ±0.053 


Ri 


20.775 ±0.027 



Table 1: Average line shape parameters from the results of the four LEP 
experiments. 



The photon contribution is taken from QED 

^ = ^Q e 2 Q?Nf (3) 
oS 

where Np is the number of colours for quarks and 1 for leptons. The inter- 
ference term contains combinations of the couplings of the fermions to the 
Z and their charge, which cannot be expressed in terms of the Z parameters 
used in equation | and are included in the Jf parameter 

4vra 2 s - Ml 

a ' mt ~ ~^s (s-M 2 ) 2 + ( S r 2 /M 2 ) 2 [ ) 

Since the interference term is expected to be small it is usually set to the 
Standard Model value introducing a small model dependence in the fit. 
The effect of initial state radiation is taken into account in the fitting pro- 
cedure through its convolution with the Born cross section. The fitted Z 
lineshape parameters are given in Table [3]. These are the physical param- 
eters, i.e they include all electroweak and strong radiative corrections and 
QED final state corrections: 



G F M| 
W2 



?A.f 



+ gv f 1 + ^QED 1 + 5qcd (5) 



The QCD correction Sqcd is zero for leptons and is to a first approximation 
proportional to (1 ±a s /7r) for hadrons. The coupling constants gA,f and gv,f 
absorb the electroweak corrections through Ap and sin 2 O 1 ^: 



A ,f = J I + Ap x I 3 g v ,f = J I + Ap x (I 3 - 2Q f sin 2 ) (6) 
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The dependence of Ap and of sin 2 O 1 ^ on M t is quadratic, while the depen- 
dence on the Higgs mass is only logarithmic giving less sensitivity to Mr. 
These radiative corrections are detailed in ||. 

The measured values of the Z parameters are confronted with the Standard 
Model prediction in order to extract the radiative corrections and therefore 
to determine a s , Mh and M t . This fit is performed using the latest version 
of the programs ZFITTER and TOPAZO || which include new calculations 
of radiative corrections: 



non factorisable QCD/EW corrections |T(| resulting in an increase of 



the fitted value of a s of +0.001 with respect to old versions. 

two loop irreducible EW corrections 0(a 2 M^/M^ v ) leading to a 
decrease in the Higgs mass of about 30 GeV/c 2 . 



• four loop QCD corrections in the (3 function ||12|| . 
Results will be given in the next Sections. 



3 Extraction of a s 

The most sensitive parameters to a s are T z , <r° ad and R^. They all depend on 
a s through Thad- Their dependence with a s is given in a first approximation 
by: 

Re = if K v 1 + t) (7) 
r z = r had + r e + r M + r r + r inv « f 1 + 0.7—) (8) 



7T 



r c rhad /. - . «i 

r had ^ T=a - - \ " 



oc (1-0.4^) (9) 

Table ^| gives the relative experimental precision on these parameters and 
the derived uncertainty on ct s when determined using only the corresponding 
parameter. R^ and Tz both allow to determine a s with the same experi- 
mental precision! 2 ] but Tz varies rapidly with Mh introducing an additional 



2 The error on cr{J ad should decrease from 1.3 x 10 -3 to about 0.9 x 10 -3 with the 
decrease of the theoretical error on the luminosity but er° ad will still be less powerful 
than Rf to determine a s . 
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Ax/x 


Aa s 


r z 

°had 


1.3 10" 3 
1.0 10~ 3 
1.3 10" 3 


0.004 
0.004 
0.010 



Table 2: Relative experimental precision on the parameters used in a s deter- 
mination and induced uncertainty on a s . 

uncertainty if no constraint on Mh is used whereas the Mh dependence of R^ 
almost cancels in the ratio of the hadronic to leptonic widths. Therefore R^ is 
the most sensitive to a s and allows to determine a s without any assumption 
on the Higgs mass. 

One should note that R^, Tz and <r° ad all depend on Tb through i\ a d and 
that R^ = r b /r had is fixed to its Standard Model value in order to extract 
a s . However, since R° is very sensitive to new physics through the Zbb ver- 
tex corrections its value may deviate from Standard Model expectation. The 
experimental determination of R° is shown in Figure |2| and is in agreement 
within 1.8cr with the SM prediction for the direct top mass measurement 
M t = (174.1 ± 5.4) GeV/c 2 [|l3|. This value of M t is used in the following. 

The measurement of R^ 

As already discussed in Section |2|, the uncertainty on R^ is dominated by 
the dilepton statistical and systematic experimental errors. The preliminary 
measurements of the 4 LEP experiments []5| presented in Jerusalem 97 are 
shown in Figure || The statistical uncertainty is AR^/R^stat) ~ (0.05(qq)© 
0.15(1 + 1~))% for each LEP experiment. The systematic uncertainty from the 
hadronic channel is of the order of 0.08% and an effort is being made to 
reduce the systematic error from the leptonic channel. 
A new measurement of R^ using ALEPH data with a global dilepton selec- 
tion allows to reduce the systematic uncertainty from the leptonic channel to 
0.08%. The dilepton selection has to be inclusive in order to be compared to 
the theoretical prediction. Therefore this selection includes also four fermion 
final state events of the type £ + £~V where i is a lepton and V is a low in- 
variant mass pair of fermions. 

First, dileptons are selected with an efficiency of 99.2% inside the detector 
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r ir 

4 K' A k, 



ALEPH mult 

1992-95 

DELPHI mult 

1994-95 

L3 mult 

1994-95 

L3 shape 
1991 

OPAL mult 

1992- 94 

SLD vtx mass 

1993- 97 



LEP leptons 



LEP+SLC 

250 




0.2159 ±0.0009 ±0.0011 
0.2166 ± 0.0008 ± 0.0009 
0.2179 ± 0.0015 ± 0.0026 
0.2223 ± 0.0030 ± 0.0064 
0.2178 ± 0.0014 ± 0.0017 
0.2158 ± 0.0017 ± 0.0014 

0.2227 ± 0.0020 ± 0.0025 



0.21736 ± 0.00086 

corrected for Yexch:jii^r 



0.215 0.22 0.225 

r^,, for ryr tad = 0.172 



Figure 2: Experimental determination ofR® at LEP and SLC. Also shown is 
the Standard Model expectation as a function of M t and the direct top mass 
measurement M t = (174.1 ± 5.4) GeV/c 2 pl 



acceptance and background from 77, qq and cosmic events is reduced to 
the level of 0.2%. All the systematic uncertainties are at the level of few 
10 -4 . Then the lepton flavour separation is performed inside the dilepton 
sample so that the systematic uncertainties are anti-correlated between 2 
lepton species and that no additional uncertainty is introduced on R^. This 
separation is needed in order to identify e + e~ — > e + e~ events for which the 
t-channel contribution has to be subtracted. This subtraction is performed 



using theoretical calculation described in ||14|| . The theoretical error assigned 
to the subtraction is 0.08% of the e + e~ — > e + e~ s-channel cross section |i~5|. 
The value of R^ obtained with this new selection is 

R £ = 20.732 ± 0.038 (ALEPH) (10) 

The relative error on R^ is reduced from 2.4 x 10~ 3 (1997 value) to 1.9 x 10 -3 
(this measurement). 
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Ratio R, = r had /r, 



> 

CD 

o 

X 

E 



10'- 



ALEPH •- 

DELPHI » 

L3 

OPAL 

LEP 



10' 




20.7 



20.754 ±0.049 
20.759 ± 0.063 
20.788 ± 0.066 
20.834 ± 0.056 

20.775 ±0.027 

X 2 /DoF: 1.4/3 



I1/a= 128.896 ±0.090 
la s = 0.118 + 0.003 
]m t = 174.1 ±5.4 GeV 



20.8 20.9 

R, 



21.0 



Figure 3: Determination of at LEP. The Standard Model prediction as 
a function of Mh is also shown. The width of the Standard Model band 
corresponds to the uncertainties on a s , M t and a(M|). The total width of 
the band is the linear sum of these effects. 



Determination of a s with R^ 

The dependence of R^ on a s is parametrised with the latest version of ZFITTER|5| : 
Rt = 19.934(l + 1-045^) + 0-94^) 2 - 15(^) 3 ) (11) 

In this parametrisation Mh = 300 GeV/c 2 and M t = 174.1 GeV/c 2 are used. 
The small dependence with the Higgs and the top masses is also parametrised 
with ZFITTER 

OX—ID) <*> 

The values of R^ obtained with this parametrisation agree with ZFITTER 
prediction at the level of few 10 -5 for any value of a s in [0.100,0.130], M t 
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in [150,200] GeV/c 2 and M H in [60,1000] GeV/c 2 . Figure | shows the new 
ALEPH determination of and the Standard Model expectation as a func- 
tion of a s . The variations arising from Mh and M t are also shown. The value 



0.725 



0.12 



0.115 

«. 

0.11 



0.105 
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20.6 20.65 20.7 20.75 20.8 20.85 



Figure 4: ALEPH determination of R? (hatched band) and Standard Model 
expectation as a function of a s . The full lines give the SM prediction for 
M t = 174.1 GeV/c 2 and M H =70, 300, and 1000 GeV/c 2 . The doted lines 
correspond to a variation ofM t of ±5 GeV/c 2 for Mh fixed to 300 GeV/c 2 . 

obtained for a s is 

a s (M 2 z ) = 0.119 ± 0.006 cxp ± 0.002 QCD ± 0.002 Mh (ALEPH) (13) 

where the first error comes from the experimental uncertainty on R^, the 
second error covers missing higher-order corrections and uncertainties in the 
interplay of electroweak and QCD corrections |16[ and the last error is ob- 
tained by varying the Higgs mass from 60 to 1000 GeV/c 2 . 
The preliminary LEP combined value or (Jerusalem 97, see Ref. 0) is 
shown in Figure |3| and leads to 

a s (M 2 z ) = 0.125 ± 0.004 cxp ± 0.002 QCD ± 0.002 Mh (14) 
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In an overall fit using all the information from LEP1 and SLC data (cross 
sections and asymmetries) the value of a s and of Mh are simultaneously 
constrained (see Section [|), the fitted value of a s is 

a s (M z ) = 0.120 ± 0.003 cxp ± 0.002 QCD (15) 

The experimental error is reduced from 0.004 (in equation |TJ]) to 0.003 be- 
cause T z and cx° ad are used in this last fit, and the 0.002 error from the Higgs 
mass disappears since Mh is also constrained. This value of a s is lower than 
in the fit to alone (equation [14]) because T z prefers lower values of a s and 
the low fitted value of M H (66 GeV/c 2 ) brings a s down by ~0.002 (in the fit 
to Re the value M H = 300 GeV/c 2 was used). This value is in good agree- 
ment with the world average a s (M|) = 0.118 ± 0.003 JT^] and of comparable 
precision. 



4 Constraint on the SM Higgs mass 

The determination of sin 2 O 1 ^ 

At the Z resonance, the most sensitive parameter to the Higgs mass is 
sin 2 which is determined through the measurement of the asymmetries 

• the lepton, b and c quarks Forward-Backward asymmetries Ap B = |«4e«4f 
determined at LEP1 from e + e~ — > ff angular distributions 

<!<7 oc 1 + cos 2 6> + |a fb cos6» (16) 



dcos6* 3 

• The t polarisation (LEP1) 

A T (1 + cos 2 6) + 2Acos0 
VAC0S9) = ~ 1 + cos 2 ^ + 2A c A T cos9 

• the Left-Right asymmetries measured at SLC with polarised electron 
beam 

A L r = A c (18) 

where At = 2(gv,f/gA,f)/(l + (gv,f/gA,f) 2 )- All these quantities are expressed 
in terms of the effective mixing angles of leptons, sin 2 6^*. The derived val- 
ues of sin 2 O 1 ^ 1 from the different asymmetry measurements are compared 
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in Figure [|. SLD and LEP average differ by 2cr and SLD data prefer lower 
values of the Higgs mass. Since all LEP data have been analysed new mea- 
surement can only come from SLC with increased statistics. 



The global fit 



A, 
\ 

<Q fb > 



Average(LEP) 
A |r (SLD) 
Average(LEP+SLD) M 



10 ° t| 



> 
o 

O 



E 10 2 



0.23 



Preliminary 

0.2310210.00056 
0.23241 ± 0.00080 
0.2319310.00090 
0.23211 ±0.00039 
0.231610.0010 
0.2321 ±0.0010 
0.2318510.00026 

X 2 /d.0.f.:3.3/5 

0.23084 ± 0.00035 
0.2314910.00021 

X 2 /d.0.f.:8.6/6 



m1/<x= 128.896 ±0.090 
raa s =0.118±0.003 
tfflm f 174.1 ±5.4 G eV 



0.232 0.234 

2 lept 



0.236 



sin 6 



eff 



< 



Figure 5: Determinations of sin 2 9^ 
from the asymmetries. The predic- 
tion from the Standard Model as a 
function of Mh is also shown. The 
width of the Standard Model band is 
due to the uncertainties in a(M|), 
a B (M|) and M t . 
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Figure 6: Ax 2 



10 

m H [GeV] 



V — V 



10 



vs. M H 



curve. The line is the result of the fit 
using all data (last column of Table 
Qj. The vertical band shows the 95% 
CL exclusion limit on Mh from direct 



searches f$DJ. 



In the global fit the electroweak measurements are compared to Standard 
Model predictions in order to determine simultaneously the values of a s (M|), 
M t and M H . The parameters G F = (1.16639 ± 0.00002) x 10~ 5 GeV~ 2 from 
muon decay and «(M|) 1 = 128.896 ± 0.090 from (JTSj are also used as input 
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in the fit. This global fit is performed with three data sets, results are given 
in Table || 

1. A first fit to LEP data alone including the measurement of the W mass 
at LEP 2 is made (first column of Table §). This fit shows that LEP 
data prefer a light top and a light Higgs. 

2. Then LEP1 and SLD data and the measurement of sin 2 #w from uN 
experiments are used to determine top quark and W masses indirectly 
(second column of Table |||). The good agreement with the direct 
measurements M t = 174.1 ± 5.4 GeV/c 2 |13| and M w = 80.430 ± 
0.080 GeV/c 2 |]| provides a test of the Standard Model. 

3. All electroweak measurements including the direct M t and Mw mea- 
surement are used to obtain the best constraint on Mh (last column of 
Table |3|) . Figure [I] shows the result of this fit using the latest version 
of ZFITTER and TOPAZ0. Also shown is the curve obtained when 
new higher-order corrections are neglected ('no (9(g 4 M 2 /M w )' curve). 
These corrections lead to a decrease of the Higgs mass of ~ 30 GeV/c 2 . 
The upper limit on Mh is 215 GeV/c 2 at 95% CL. The uncertainty on 
a(M|) causes an error of 0.2 on Log(Mn). The fit is also performed 
with the new evaluation of a(M|)P0|. With this value the fitted error 
on Log(M H ) is reduced by 30% (from 0.33 to 0.25). 

Since Mh is mainly constrained by the measurement of sin 2 ^^*, new data 
from SLC will improve this constraint. Moreover, radiative corrections are 
shared between the top quark and the Higgs, therefore the reduction of the 
error on the direct top mass measurement will improve the constraint on the 
Higgs mass. The W mass is also sensitive to the Higgs mass and its mea- 
surement with an error of 30 MeV at LEP2 will also improve the constraint 
on Mh- 



5 Conclusion 

The large sample of data accumulated at LEP1 allows the determination of 
a s through radiative corrections, a s = 0.120 ± 0.003 CX p ± 0.002qcd with a 
precision comparable to the world average (in which LEP1 data have not 
been included). Since all LEP1 data have been analysed, this determination 
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LEP including 
LEP2 M w 


all data except 
M t and M w 


all data 


M t [GeV] 
M H [GeV] 
Log(M H /GeV) 
«s(Ml) 


157^ 
56^ 

i 7 cr+0.45 
L - ' °-0.35 

0.122 ±0.003 


16111 
33lf 7 

i cq+0.37 
l.dO_ 29 

0.121 ±0.003 


171.1 ± 5.1 

66t 3 7 9 

0.120 ±0.003 


X 2 /d.o.f. 


6/9 


14/12 


17/15 


sin 2 fl$ 

sin 2 6\y 
M w [GeV] 


0.23183 ±0.00025 
0.2245 ± 0.0008 
80.302 ± 0.040 


0.23145 ±0.00021 
0.2235 ±0.0008 
80.351 ± 0.040 


0.23146 ±0.00022 
0.2230 ±0.0005 
80.380 ±0.027 



Table 3: Results of the electroweak fits. See text for details. 

may only slightly improve through the reduced error on cr^ ad [0. 
The global fit gives for the Higgs mass 

M H = 66l^ GeV/c 2 (19) 

There is still room for improvement in this constraint through better top and 
W mass direct measurements, new data from SLC and the reduction of the 
error on a(M|). 

A low value of the Higgs mass is preferred by the Z resonance data: it is still 
possible to find the Higgs at LEP2 ! 
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